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ABSTRACT 

We present the discovery of an outflowing ionized wind in the Seyfert 1 Galaxy Mrk 335. 
Despite having been extensively observed by most of the largest X-ray observatories in the last 
decade, this bright source was not known to host warm absorber gas until recent XMM-Newton 
observations in combination with a long-term Swift monitoring program have shown extreme flux 
and spectral variability. High resolution spectra obtained by the XMM-Newton RGS detector 
reveal that the wind consists of three distinct ionization components, all outflowing at a velocity 
of ^5000 km/s. This wind is clearly revealed when the source is observed at an intermediate flux 
state (2-5x 10 -12 ergs cm~ 2 s _1 ). The analysis of multi-epoch RGS spectra allowed us to compare 
the absorber properties at three very different flux states of the source. No correlation between 
the warm absorber variability and the X-ray flux has been determined. The two higher ionization 
components of the gas (log£~2.3 and 3.3) may be consistent with photoionization equilibrium, 
but we can exclude this for the only ionization component that is consistently present in all flux 
states (log£~1.8). We have included archival, non-simultaneous UV data from HST (FOS, STIS, 
COS) with the aim of searching for any signature of absorption in this source that so far was 
known for being absorption-free in the UV band. In the COS spectra obtained a few months 
after the X-ray observations we found broad absorption in CIV lines intrinsic to the AGN and 
blueshifted by a velocity roughly comparable to the X-ray outflow. The global behavior of the 
gas in both bands can be explained by variation of the covering factor and/or column density, 
possibly due to transverse motion of absorbing clouds moving out of the line of sight at Broad 
Line Region scale. 

Subject headings: Active Galaxies: general — Active Galaxies: Mrk 335 

Introduction of many AGN both at UV and X-ray wavelengths 

(see Crenshaw ct al. 2003 for a review). 



Outflowing photoionized gas is a well-established 
feature of over half of all Seyfert Galaxies. The 



There arc several reasons for keeping our in- 
terest high in the study of warm absorbers. The 
presence of this gas is revealed as a series of , , , , , , , , , , 

f, , , . . . . observed velocity shift (almost always to the blue) 

blueshifted absorption signatures m the spectra 
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provides evidence that material is traveling out- 
ward from the central region of AGN. If this ma- 
terial eventually leaves the AGN, then outflows 
might carry significant mass out of the AGN and, 
as a consequence, give a substantial contribution 
to the chemical enrichment of the intergalactic 
medium (IGM) (Furlanetto et al. 2001, Cavaliere 
et al. 2002, Germain et al. 2009, Hopkins et al. 
2010, Barai et al. 2011). The outflow can also 
impact the development of the host galaxy itself. 
If it is as strong as 0.5-5% of the Eddington lu- 
minosity of the AGN, then the feedback from the 
AGN can regulate the growth of the galaxy and 
the growth of the central black hole as well (Di 
Matteo et al. 2005, Hopkins et al. 2010). 

The fact that X-ray absorbers are always found 
in sources with UV absorption (Crenshaw et al. 
2003, Kriss 2006) indicates that there is some in- 
terplay between the two phenomena and this in 
turn implies that photoionized gas must be dis- 
tributed with high coverage in the AGN. 

Several open questions are still awaiting solu- 
tion. In many cases, high quality observations of 
bright sources have allowed us to relate kinemat- 
ically blue-shifted X-ray and UV absorption lines 
(e.g. Kaspi et al. 2002, Kaastra et al. 2002) lead- 
ing to the idea of an outflowing wind that spans 
a wide range of ionization states with a common 
velocity pattern. In contrast, other cases have 
shown no obvious relation between UV and X- 
ray absorption components (Kriss et al. 2011). 
The ionization structure of the wind was some- 
times found to be a continuous distribution, as 
proposed for NGC 5548 (Steenbrugge et al. 2005) 
or made by discrete, separate ionization compo- 
nents, sometimes in pressure equilibrium (Kron- 
gold et al. 2003), and other times not (Detmers et 
al. 2011). The location of the photoionized gas is 
also a subject of debate. The Narrow Line Region 
on kpc scales (Kinkhabwala et al. 2002), the inner 
edge of the molecular AGN torus on pc scale (Kro- 
lik & Kriss 2001) and the accretion disk on sub-pc 
distances (Elvis 2000, Krongold et al. 2007) have 
all been put forward as possible regions of origin 
for these winds. 

One sound piece of evidence that so far has been 
fully supported observationally is that, in the vast 
majority of sources -practically 100% - intrinsic 
absorption is present in both X-ray and UV bands 
(Crenshaw et al. 2003, Kriss 2006). 



The latest XMM-Newton and Hubble Space 
Telescope (HST) observations of the Narrow Line 
Seyfert 1 Mrk 335 (z = 0.025785, Huchra et al. 
1999) presented in this paper resolve a potential 
challenge to this paradigm. 

Mrk 335 has a long X-ray history. It has been 
observed by all the largest facilities throughout 
the years, with the exception of Chandra. The X- 
ray spectrum has always shown a standard AGN- 
shapc: mildly steep power law in the 2-10 keV 
band, a strong soft X-ray excess and a complex 
and prominent Fe K line (see Larsson et al. 2007, 
Longinotti et al. 2007, O'Neill et al. 2007). 

Intrinsic soft X-ray absorption was never de- 
tected with high confidence. In the ASCA survey 
of warm absorbers presented by Reynolds (1997), 
this source was not among the list of the ones pre- 
senting the characteristic absorption edges that re- 
vealed the presence of ionized gas along the line of 
sight. Likewise, Mrk 335 has never shown any UV 
absorption either (Zheng et al. 1995, Crenshaw et 
al. 1999, Dunn et al. 2007). 

In 2007, Mrk 335 underwent an abrupt decrease 
in flux (Grupe et al. 2007, 2008) and it gave us the 
opportunity to observe and study a previously un- 
known emission-line component that revealed the 
presence of photoionized gas at Broad Line Region 
scales (Longinotti et al. 2008). The broadband 
(CCD) spectrum of this peculiar low flux state of 
Mrk 335 was explained by Grupe and collabora- 
tors by an intervening partial covering absorber, 
although a competing and alternative scenario in 
which the low state is produced by disk reflection 
could not be rejected. Noticeably, soft X-ray ion- 
ized absorption was not observed in the high res- 
olution data of any flux state (O'Neill et al. 2007, 
Longinotti et al. 2008). 

In 2009 Mrk 335 was observed for 200 ks by 
XMM-Newton as part of a trigger program based 
on Swift long term monitoring (Grupe et al. 2012). 
The source was caught for the first time at a flux 
state intermediate between the deep minimum of 
2007 and the "standard" bright flux in which it 
was always observed. These data present very 
different characteristics with respect to previous 
XMM-Newton observations. In particular they 
clearly show intrinsic soft-ray absorption in the 
RGS band. This finding seems to question the 
one-to-one correspondence reported between X- 
ray and UV absorbers. 
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Table 1: XMM- Newton Observation log. 1) Gon- 

doin et al. 2002; 2) Longinotti ct al. 2007; 3) O' Neill et 
al. 2007; 4) Grape ct al. 2008; 5) Longinotti ct al. 2008; 
6) Grape ct al. 2012; 7) Gallo et al. 2013 



OBSID 


Date 


Exp 


Flux state 


Rcf 






(ks) 






0101040101 


2000-12-25 


37 


High 


1,2 


0306870101 


2006-01-03 


133 


High 


3 


0510010701 


2007-06-11 


23 


Low 


4,5 


0600540601 


2009-06-11 


132 


Mid 


6,7 


0600540501 


2009-06-13 


82 


Mid 


6,7 



The present paper presents the first high- 
resolution study of the warm absorber in Mrk 
335 and it is based on the comparison of all the 
data obtained up to the year 2009 by the RGS 
spectrometer onboard XMM-Newton. 

The X-ray analysis presented here is com- 
plemented by the inclusion of recent and non- 
simultaneous HST-COS spectra that were ob- 
tained a few months after the X-ray observations. 

2. Observations and data reduction 
2.1. XMM-Newton data 

Mrk 335 was observed 5 times by XMM-Newton 
(see Table [T] for the log of observations). All the 
data sets were reduced with SAS 11.0.(0. For a de- 
tailed description of the various data sets, we defer 
the reader to the numerous prior publications on 
this source that are listed in Table [T] 

The focus of the present paper is the analysis of 
the warm absorber features in the high-resolution 
spectra obtained by the Reflection Grating Spec- 
trometer (RGS, den Herder et al. 2001), therefore 
we do not include CCD data from the EPIC MOS 
cameras, which are extensively described in pre- 
vious publications. However, the EPIC pn data 
will be shown and used for measuring fluxes and 
constraining models, therefore we briefly describe 
their data reduction in the following. 

For all 5 data sets, the EPIC pn camera 
(Struder et al. 2001) was used as prime instru- 
ment, and it was operated in Full Frame mode 



http: / /xmm.esac.esa.int/sas/current/documcntation/sas_concise.; 
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Fig. 1.— Multi-epoch EPIC pn spectra of Mrk 
335. The effective area has been accounted for, 
therefore the difference between the various obser- 
vations has to be attributed to intrinsic changes in 
the spectral shape. 

(0101040101, 0600540501, 0600540601), Small 
window mode (0306870101) and Large window 
mode (0510010701). Data reduction was per- 
formed by using the standard tasks epproc and 
rgsproc, for pn and RGS data respectively. 

Screening of high particle background was ap- 
plied to the pn data so as to maximise the signal- 
to-noise ratio according to the procedure described 
in Piconcelli et al. (2004) . The spectral extraction 
regions in the pn camera is a circle of 40 arcsec of 
radius for source and background spectra, in all 
observations. 

2.1.1. The EPIC pn data 

The CCD data of the XMM-Newton campaign 
are the subject of two more publications presented 
by our group (Grupe et al. 2012 and Gallo et 
al. 2013). The former paper reports on the long 
term variability of Mrk 335, which has been moni- 
tored by the Swift X-ray telescope since 2007, and 
the application of the partial covering model to 
explain the variation observed by XMM-Newton. 
The latter paper concentrates instead on the anal- 
ysis of the disk reflection component in the XMM- 
Newton 2009 data. 

This complementary work provides the first 
detailed description of the line-of-sight gas in 
Mrk 335 as seen by a high resolution spectrome- 



Table 2: X-ray fluxes (in units of 
iCP 11 erg cm~ 2 s _1 ) measured by fitting the 
EPIC pn data with a broken power law model 
with break energy fixed at 2 keV. 



OBSID 


Flux (0.3-2 kcV) 


Flux (2-10 kcV) 


0101040101 


2.60± 0.01 


1.37± 0.01 


0306870101 


3.32± 0.01 


1.77± 0.01 


0510010701 


0.20± 0.01 


0.34± 0.01 


0600540601 


0.36± 0.01 


0.48± 0.04 


0600540501 


0.52± 0.02 


0.51± 0.05 



ter, and, as to the X-ray band, it is entirely based 
on the RGS data. However, consistency between 
spectral fits performed on the RGS and the EPIC 
pn data is explored for completeness. 

To generally illustrate the dramatic change in 
the spectral shape of Mrk 335, we have included all 
the EPIC-pn spectra of Mrk 335 in Figure [U The 
0.3-10 keV figure shows very effectively that the 
spectral changes are concentrated below 4-5 keV, 
and that the difference at higher energy is merely 
attributed to a normalization effect. Continuum 
fluxes extracted from these spectra are listed in 
Tabled 

2.2. HST Ultraviolet Spectra 

A few months after the XMM-Newton obser- 
vations of Mrk 335 in 2009, the Cosmic Origins 
Spectrograph (COS) team observed Mrk 335 as 
part of their program to probe warm and hot gas 
in and near the Milky Way using AGN as back- 
ground sources (Proposal ID 11524, James Green, 
PI). Green et al. (2012) describe the key char- 
acteristics of the design and performance of the 
COS instrument on the Hubble Space Telescope 
(HST). We have retrieved these spectra from the 
HST archive to see if any trace of the unexpected 
X-ray absorption we see in the RGS spectra has a 
UV counterpart. Table [3] summarizes the observa- 
tional details. 

The first observation on 2009 October 31 used 
the Primary Science Aperture (PSA) and grating 
G160M with four different central wavelength set- 
tings. All four observations used the same focal 
plane position for the detector, FPPOS=3.11 four 
observations used the same focal plane position 
for the detector, FPPOS=3. The second obser- 



vation on 12 February 2010 used gratings G130M 
and G160M, again with four central wavelength 
settings for each grating and FPPOS=3. The four 
different central wavelength settings provide full 
coverage across the wavelength gap between seg- 
ments A and B of the FUV detector. They also 
place the grid- wire shadows and other detector ar- 
tifacts at independent places along the spectrum 
so that they can be more easily removed from the 
data without gaps in wavelength coverage. Wc 
used vl2.17.6 of the COS calibration pipeline to 
process and combine the data. Our processing 
included a 1-dimcnsional flat-field correction that 
eliminates the grid-wire shadows and other detec- 
tor artifacts. The combined G160M spectra cover 
the wavelength range 1400-1770 A; the G130M 
spectrum covers 1135-1440 A. The G160M spec- 
tra each have a median signal-to-noise ratio (S/N) 
of ~ 17 per resolution element. For G130M, the 
median S/N is ~ 33 per resolution clement. 

Although the COS line-spread function (LSF) 
has broad wings that can fill in narrow absorp- 
tion features (Ghavamian et al. 2009), we do not 
detect any intrinsic lines in Mrk 335 that would 
necessitate deconvolving the effects of the LSF as 
in Kriss et al. (2011). Ghavamian et al. (2009) 
note that the broad LSF has negligible impact 
on spectral features with Doppler widths exceed- 
ing ~ 100 kms~ . All spectral features of inter- 
est directly related to Mrk 335 have much greater 
widths. 

After processing through the pipeline, we mea- 
sured the locations of strong interstellar lines in 
the spectra to determine a zero-point correction 
for the wavelength scale using the H I velocity of 
Vlsr = —11 km s _1 along the Mrk 335 sight- 
line (Murphy et al. 1996). Overall uncertainties 
in the COS wavelength scale limit our knowledge 
of the absolute wavelength scale to siml5 kms -1 . 
Figure [2] shows the calibrated COS spectrum of 
Mrk 335, with the G130M observation from Febru- 
ary 2010 merged with the G160M observation of 
October 2009. 

To place these most recent HST spectra in an 
historical context, we also retrieved archival spec- 
tra obtained using the Faint Object Spectrograph 
(FOS) in 1994, and more recent Space Telescope 
Imaging Spectrograph (STIS) spectra obtained in 
2004. The relevant data sets are also listed in Ta- 
ble [3] For these data we performed no special 
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Table 3: HST Observations of Mrk 335 



Data Set Name 


Inst runic nt 


Grating/ Tilt 


Date 
(GMT) 


Start Time 


Exposure Time 


Ib4n05010 


cos 


G160M/1589 


2009-1 0-31 


13:40:19 


408 


Ib4q05020 


COS 


G160M/1600 


2009-10-31 


13:50:24 


408 


Ib4a05030 


COS 


G160M/1611 


2009-10-31 


14:00:29 


408 


Ib4n05040 

1 Ut:U V 7 7 V 7 1 V 7 


COS 


G160M/1623 


2009-10-31 


14:10:34 


409 


ih4nnfim n 

1 U4:U VJUVJ X VJ 


COS 


G1 finM /I c >89 

vjr 1UU1V1 / 100-7 


201 0-02-08 

ZiVJ-LVJ VJZi VJO 


07-34-24 


302 

tjVJZi 


Ib4n06020 


COS 


VJX J- VJ VJ IV J. / X VJ VJ VJ 


2010-02-08 

ZJVJ-LVJ VJZi vju 


07-43-02 

VJ 1 . t:(J . VJZj 


302 


ib4n06030 

1 L»^tV_^VJVJVJOVJ 


COS 


GlfiniVT/lfil 1 

VJX 1UU XYJ- / _L VJ _L _L 


2010-02-08 

ZjVJJ-VJ VJ Lj VJO 


07:51:40 


302 


Ib4q06040 


COS 


G160M/1623 


2010-02-08 


08-00-18 


301 


ih4n06050 


COS 


G130M/1291 

VJT i (JVJ XV _L / _L _ !7 -L 


2010-02-08 

ZjVJJ-VJ \J iu VJCJ 


08-58-05 


608 


Ih4n06060 

1 U^tV^VJVJ VJUVJ 


COS 


O130M/1300 

vjr xovjivx / xovjvj 


2010-02-08 

ZjVJ-LVJ \J iu VJCJ 


09-1 5-01 

VJ ij . X (J . VJ X 


605 


lU^tV^VJVJVJ I VJ 


COS 


G130M/1309 

vjr x o vj iv j. / _i_ o vj 


2010-02-08 

ZjVJJ-VJ VJZj VJCJ 


09-28-52 


604 


ID^qUOUoU 




vjrlOUlVI/ lOlo 


ZUlU-UZ-Uo 


1U.O I .00 


DUO 


y29e0202t 


FOS 


H13 


1994-12-16 


05:43:31 


1390 


y29e0203t 


FOS 


H13 


1994-12-16 


07:01:06 


770 


y29e0204t 


FOS 


H19 


1994-12-16 


07:20:06 


960 


y29e0205t 


FOS 


H27 


1994-12-16 


07:40:07 


60 


y29e0206t 


FOS 


H27 


1994-12-16 


08:40:19 


420 


y2gq0301t 


FOS 


PRISM 


1994-12-16 


08:48:35 


60 


o8n505010 


STIS 


E140M 


2004-07-01 


16:24:51 


1945 


o8n505020 


STIS 


E140M 


2004-07-01 


17:44:50 


2295 


o8n505030 


STIS 


E140M 


2004-07-01 


19:20:49 


2875 


o8n505040 


STIS 


E140M 


2004-07-01 


20:56:49 


2875 
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Fig. 2 — Calibrated COS spectrum of Mrk 335. Wavelengths shortward of 1430 A are from the G130M ob- 
servation of February 2010; longer wavelengths are from the G160M observation in October 2009. Prominent 
emission features are indicated. Geocoronal emission in the center of the galactic Lya absorption trough 
is indicated with an Earth symbol. The narrow absorption features in the spectrum are either foreground 
interstellar lines or intervening intcrgalactic Lya absorbers. 
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processing other than to correct the zero point of 
the wavelength scale of the FOS spectra to align 
the ISM absorption features with the H I veloc- 
ity noted above. The STIS spectrum required no 
additional correction. Table 0] compares the con- 
tinuum fluxes for all the HST observations of Mrk 
335 that we discuss below. 

3. X-ray spectral analysis 
3.1. Strategy 

The spectral analysis of the XMM- Newton RGS 
spectra was carried out by using two indepen- 
dent photoionization codes. We applied the xabs 
warm absorber model included in the SPEX fitting 
package (Kaastra, Mewe, Nieuwenhuijzcn, 1996). 
In parallel, the spectra were analyzed following 
the same procedure by employing the photoion- 
ization code PHASE, developed by Krongold et 
al. (2003). All the results presented in this pa- 
per are extracted from the analysis carried out 
with the SPEX software. However, the warm ab- 
sorber properties derived for the 2009 spectra and 
reported in the following section were fully con- 
firmed by the second analysis. In this way we are 
confident to have checked the possible model de- 
pendency that sometimes is introduced by the use 
of different photoionization codes. 

During the fitting procedure RGS data have 
been binned by a factor of 4. Whenever other 
choices of grouping were taken (i.e. for plotting 
purpose), it will be specified. x 2 -statistics was ap- 
plied and 1-sigma error bars are quoted through- 
out the paper. 

We highlight our strategy in the following. 
Since the main purpose of the present work is 
a detailed study of the warm absorber and its 
multi-epoch behaviour, we start by analysing the 
data of 2009, where the presence of the ionized 
absorber is very evident. 

Once a complete model of the ionized gas is 
established on sound bases for the mid state, we 
extend it to the other data sets. In this way we 
will test the behaviour of the absorber at different 
epochs and its relation to the corresponding flux 
levels that Mrk 335 has shown during the XMM- 
Newton observations. 

Figure [3] displays the comparison of the RGS 
residuals of the mid and high state data fitted by 



a simple power law and Galactic absorption model. 
The absorption around 16 A, which is obvious in 
the mid state spectrum observed in 2009, is un- 
doubtedly less pronounced in the high state data 
of 2006. Analogously, the emission line around 
19 A, is more evident in the mid state data, due 
to the lower continuum flux with respect to the 
2006 spectrum. After this first-order comparison, 
we proceed to model the absorption observed in 
the mid state spectra. 

3.2. The ionized gas in the 2009 mid state 

We start by fitting the continuum in the RGS 
band (7-30 A) with a power law and a black body 
component to mimic the soft excess. The Galactic 
column density along the line of sight to Mrk 335 
is also included and fixed to 4xl0 20 cm -2 (Dickey 
& Lockman 1990). 

To this purely phenomenological model of the 
continuum we added the effect of absorption from 
ionized material. This effect is modeled by the 
xabs warm absorber code that reproduces the 
transmission of the nuclear continuum through a 
slab of ionized gas in photoionization equilibrium. 
The free parameters in our fits are the column den- 
sity of the gas, the ionization parameter, the out- 
flow velocity and the root mean square velocity 
width of the absorption lines. Solar System abun- 
dances from Lodders et al. (2009) are assumed. 
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Mid stote 0600540601 11/06/2009 

Mid state 0600540501 13/06/2009 




0.100 
Energy (keV) 



Fig. 4. — Observed spectral energy distribution 
of Mrk 335 in the 2009 mid state from XMM- 
Newton. 

Initially, the ionized absorber is assumed to be 
totally covering the source of radiation. Subse- 
quently, in the fitting procedure we will test for 
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Table 4: Continuum Fluxes for HST Observations of Mrk 335 



Observatory 


Date 


F(1500 A) 






(10" 14 ergs enr 2 s~ x A X ) 


HST/FOS 


1994-12-16 


6.0 


HST/STIS 


2004-07-01 


3.6 


HST/COS 


2009-10-31 


3.1 


HST/COS 


2010-02-08 


3.1 



the possibility of a more patchy geometry of the 
gas. 

The absorption spectrum is calculated assum- 
ing the photoionization balance produced by the 
incoming ionizing radiation. 

Figure |4] shows the spectral energy distribu- 
tion (SED) of Mrk 335 including also the data 
from the Optical Monitor (OM) onboard XMM- 
Newton. The OM observed Mrk 335 in 4 filters 
centered at 2182, 2341, 2946 and 3481 A, simul- 
taneously to the X-ray instruments. This allowed 
us to add the information on the ultraviolet part 
of the SED for calculating the ionization balance 
in our warm absorber models. 

Figure Q] and Table [2] show a few percent differ- 
ence in the soft X-ray flux of the two observations 
of 2009 (see also Grupe et al. 2012 for the light 
curve analysis). A visual inspection of the data 
also reveals little difference in the spectral shape 
of the two RGS spectra. We therefore decided to 
start by carrying out the analysis of the two data 
sets separately and check consistency in the spec- 
tral parameters during the fitting process. All the 
fit statistics reported below refer to the data set 
extracted from OBSID 0600540501. 

A continuum model without any intrinsic 
absorption consisting of Nh Gal * (power law 
+ black body) produces an unsatisfactory fit 
(x 2 /d.o.f. = 1516/1052). 

The addition of one warm absorber compo- 
nent to this baseline model is highly significant 
(x 2 /d-O.f.=1435/1048) but still not sufficient to 
provide a good fit to the spectra. We added a 
second component and found an improved best fit 
model with xVd- - 1 ^ 1334 / 1044 - When the pres- 
ence of a third warm absorber was tested the final 
fit statistics of x 2 /d.o.f.=1307/1040 was reached. 
The data fitted by this model are shown in Fig- 
ure El The presence of an additional absorption 
component that might be postulated after looking 



at the residuals, was tested and found not signifi- 
cant. 

The best fit parameters for the three warm ab- 
sorbers are summarised in Table [SJ The separate 
contribution of each absorption component to our 
global model is illustrated in Figure|6] The 3 mod- 
els have been renormalized so as to highlight the 
dominant features in each of them. 

3.2.1. Consistency between the two mid state 
spectra of 2009 

To carefully check the consistency between 
the two separate XMM-Newton observations 
we applied the warm absorber model described 
above to the second mid state data set (OBSID 
0600540601). This spectrum is shown in the top 
panel of Figure [3] In this fit, all the parame- 
ters are left free to vary, including those of the 
baseline continuum model. The detailed results 
are reported in the second row of Table [5l We 
concluded that the fit of the two RGS data sets 
yield parameters fully consistent within the errors. 
Furthermore, we note that the soft X-ray flux vari- 
ations between the two observations of 2009 (the 
2-10 keV flux being practically constant, see Ta- 
ble [2]) , does not present variability higher than a 
30% level, which is approximately the accuracy 
in the errors of the detected column densities and 
ionization parameters (Tabic [5]). Therefore, warm 
absorber variations induced by the observed flux 
variability would not be detectable in these data 
sets. 

Nonetheless, we need to consider the results of 
the variability behaviour of Mrk 335 reported by 
Grupe et al. (2012). These authors have split the 
data of 2009 according to the pn counting rate and 
they have obtained two spectra corresponding to 
"faint" and "bright" portions within the 2009 mid 
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Table 5: Best fit parameters of the 3 warm absorber model in the mid states of 2009. The covering factor of 
the absorbers is fixed to 1. The emission lines reported in Section [3.31 are included in these fits. 



Data set 


continuum 
(T in kcV) 


Wa Phase 


(erg cm s — 1 ) 


Column density 
(cm" 2 ) 


(km s~ 1 ) 


Vbroad 

(km s" 1 ) 


X 2 /d.o.f. 


0600^40^01 


r-2 45+ ' 10 


I 


1 92 + 05 


2 72+ 1 ' 11 x 1 21 
^' ,z -1.03 x lu 


4300+ 950 


<30 


1 ^07 /I 040 




kTi,(,=0.09±0.02 


II 


2 35+ 024 


4 43+ 2 01 xlO 21 

1.30 


5450+ 1 ™ 

^ -180 


120+ 140 

-50 








III 


Q Q1 +0.09 

°- ol -0.0S 


2.13± 6 1 2 1 5 2 x lt)22 


4000+ 230 


170+ 86 




0600540601 




I 


i 7 = + 0.21 


3.27i 2 ;°Jxio 21 


60 00+ 1400 


<70 


1249/1033 




kTM,=0.11±0.07 


II 


9 I c + 0.10 


4.80l 2 ; 3 30 xl0 21 


49001 1™ 


1201™ 








III 


3 - 22 -0.3S 


2.04lf 82 xl0 22 


5200t^° 


<650 




First 80 ks 


F — 9 3O+0.30 

i — *.oo_ 5Q 


I 


1.61±0.16 


3.34±*;f° X10 21 


4700+ 850 


<20 


1226/1033 


0600540601 


kTi,i,=0.10±0.08 


II 


2 16+ 009 


4.07t 2 «xl0 21 


4900ti™ 


<70 








III 


3 26+ 17 


4.02+3 8 00 xl0 22 


52001™° 


90l™ 




combined RGS 


r=2.57±i;Ji 


I 


1 S5+ 004 


2.15t° ;f 6 xl0 21 


40001^° 


<10 


2143/1678 




kT i)i ,=0.11±0.01 


II 


9 i , +0.02 
z - 1J --0.03 


4.221°, 'JgXlO 21 


51401"" 


140±30 








III 


3 ig+0.06 
°- ly -0.09 


3.56li°|xl0 22 


5300lf 00 


120lf 




0600540601 
(RGS+ pn) 


l i -, -7 C +0.02 

r=i. 76t 002 


I 


i oQ+0.06 


2.86l^ 8 xl0 21 


47001™ 


<15 


1704/1284 




kTi,i,=0.11±0.01 


II 


n 1 c+003 
z - lo -0.02 


3.62l°; 22 xlO 21 


5100l» 


120l« 








III 


3.31±0.13 


6.22t 2 ; 2 ^xl0 22 


6200l™o 


3oi;° B 
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Fig. 3. — Comparison between residuals (in terms of a) of the mid state (top) and high state (bottom) RGS 
spectrum (OBSID 0600540601 and 0306870101) fitted only by a power law and Galactic absorption model. 
The broad absorption feature produced by the Fe UTA transitions that is very prominent in the mid state 
disappears almost completely in the high state data: residuals around 16A in the bottom panel are still 
negative but at a lower significance. 
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Fig. 6. — The plot shows the 3 (rescaled) warm ab- 
sorber components detected in the mid state spec- 
trum and reported in Table[5j From top: medium, 
high and low ionization. 



state observation^. The "faint spectrum" corre- 
sponds to the first 80 ks of obs 0600540601 and it 
is made by events with less than 3 counts per sec. 

In order to exclude that the warm absorber fit 
is driven by the brightest part of the 2009 data, 
we have applied the three warm absorbers best fit 
model to the RGS spectrum extracted from these 
80 ks. The parameters of this fit are reported in 
the third row of Table [5] In this way we test for 
possible variations of the warm absorber proper- 
ties in correlation with the source flux within the 
2009 data. The parameters derived in the "faint" 
80 ks spectrum are remarkably consistent with 
those derived when using the two integrated spec- 
tra of obs 0600540601 and 0600540501 (first two 



2 In this context "faint and bright" do not have reference to 
the multi-epoch behaviour of the source 
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Fig. 5.— The RGS spectrum of the 0600540501 mid state data fitted by the best fit model in Table [5] 
The labels of the absorption lines (grey) are blueshifted by an average outflow velocity of 5000 km/s, the 
labels for the emission lines (red) correspond instead to the laboratory values. Instrumental features can 
be distinguished by their typical "squared" shape. The relative contribution of the three warm absorbers is 
best seen in Figure El 



rows in the table). 

We conclude that the ionized gas does not 
present significant variation while XMM-Newton 
observed Mrk 335 in 2009, therefore we applied the 
best fit model to the combined data sets of 2009 
and from now on, we consider this global data set 
as "the mid state" . The results from this latest fit 
are also included in Table [SJ 

Last, we have checked the effect of the broad- 
band X-ray continuum on our warm absorber 
model. Indeed, extending the spectral band up 
to 10 keV might modify the parameters since the 
power law slope may not be correctly measured 
given the fact that the RGS band is limited to the 
soft X-rays (below 2 keV). To this purpose, we 
have applied the warm absorber model to the RGS 
and pn data simultaneously. This test was done 
for OBSID 0600540601. The resulting photon in- 
dex is indeed modified, being r=1.74+o gg, but the 
warm absorbers properties are basically consistent 
with the RGS-only based fit (see Table [5]). 



3.3. The emission lines in 2009 

When Mrk 335 was observed at its lowest flux 
state, Grupe ct al. (2008) and Longinotti ct 
al. (2008) discovered an important emission line 
component. Residuals in Figure [3] show a clear 
emission line corresponding to the position of 
OVIII Lya. 

To replicate the analysis carried out by Longinotti 
et al. (2008), we fitted the mid state data by 
adding to the three warm absorbers model a se- 
ries of emission lines with 0-width at the wave- 
lengths listed in Table HI The wavelengths were 
kept frozen to the laboratory values in order to 
allow a direct comparison to the fluxes reported 
by Longinotti et al. 2008 in their Table 2. Final 
detections are listed in Table |U Note that the 
warm absorber parameters reported in Table [5] 
and the model plotted in Figure [5] include these 
emission lines. 

Figure [7] shows the fluxes of the emission lines 
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Table 6: Emission lines detected in the mid-state 
spectra. The lines were modeled with a zero- width 
profile. A\ 2 is quoted for one free parameter. 



ID 


A 
(A) 


Flux 

X 10~ 4 phot cm -2 s _1 


Ax 2 


OVIII Lya 
OVII (r) 
OVII (i) 
OVII (f) 

NVII Lya 


18.969 
21.600 
21.801 
22.097 
24.779 


n 181 +0.050 

n 171 +0-076 
-0.060 

0.088 in™ 

n 1 ia+0.107 
U.118_ . 099 


33 
14 
12 

3 

3 



from this work (black points) and from Longinotti 
et al. (2008) (red points). From this figure, it is 
clear that there is no evidence for variability of 
the lines' intensity. It is interesting to note that 
the unusual line ratio of the OVII triplet that led 
Longinotti et al. (2008) to propose an origin in 
high density gas at BLR scale for the line emit- 
ting gas, is partially confirmed by the new obser- 
vations. The line ratios measured in the two flux 
states are consistent within 1 a and, tentatively, 
this still indicates that the emission lines arise in 
a high density gas. 

In the low state spectrum of 2007, it was possi- 
ble to measure the width of the OVIII Lya line 
(FWHM=0.14±0.05 A, corresponding to -2200 
km s _1 , see Longinotti et al. 2008). This time 
in the mid-state spectrum we could only obtain a 
loose upper limit on the width of the same line 
of 0.58 A, therefore fully consistent with the 2007 
measurement, but still too unconstrained to allow 
further speculation on the properties of the emit- 
ting gas in 2009. 

3.4. Multi-epoch variability and covering 
factor 

Once the best fitting warm absorber model is 
well established in the mid state data, we will 
take it as the reference model for tracing the be- 
haviour of the absorber along the other spectral 
states corresponding to prior epochs. Our aim is 
to investigate the relation between the ionized gas 
properties and the flux changes of the source. We 
propose to carry out this investigation by varying 
the minimum number of parameters in the best 
fit model when this is applied to the other data 
sets. This way of proceeding shall allow us to dis- 
tinguish which physical parameter is driving the 
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Fig. 7. — Comparison between the fluxes of the 
emission lines detected in the mid-state spectrum 
(black squares) and in the low-state data (red 
squares). The low state fluxes are taken from 
Longinotti et al. (2008). Data points are slightly 
offset along the x-axis for plotting purpose. 

warm absorber variability, if any. 

We shall note here the results reported in Sec- 
tion I3.2.1I XMM- Newton observed the source in 
mid state within 2 days in 2009. Despite the mod- 
erate flux change, the warm absorber was detected 
in both spectra with no appreciable change of its 
properties, but with the present data we cannot 
exclude that small variations took place on a time 
scale of days. The consistency of the model found 
when fitting the two spectra separately (see pa- 
rameters in Table [5]) can only allow us to say that 
we do not detect short term variability within the 
errors. 

Before extending our analysis to the other data 
sets, a final test is carried out on the mid-state 
spectrum. 

3.4-1- The mid state spectrum in 2009: a partial 
covering warm absorber? 

As a first consideration, we note that so far we 
have assumed that the ionized gas is distributed 
to completely cover the X-ray source. In order to 
test the possibility of a partially covering warm 
absorber, we have left the covering fraction of the 
three components free to vary in the fit of the mid 
state data. The resulting covering fraction for the 
low, medium and high ionization phases are re- 
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spectively 80, 75 and 90%, but they arc all con- 
sistent with 100%. We could measure error bars 
only for the covering factor of the low ionization 
absorber, which is constrained to 80^35%. Hence, 
the mid-state data are consistent with the pres- 
ence of a partial covering warm absorber, but this 
feature is not formally required. 

3.4-2. The high state spectrum in 2006 

The high-resolution spectrum extracted from 
this epoch at first glance does not suggest the 
presence of significant absorption (bottom panel 
in Figure [3]) , and it is well described by a sim- 
ple power law with T=2.75±0.01. We applied the 
three warm absorbers model of the 2009 mid state 
to these data assuming zero column density, i.e. 
equivalent to having no absorption. The resulting 
fit gives x 2 /d.o.f. =1376/1063. If all the absorbers 
parameters are fixed to the values reported for the 
mid state in Tabled we found, not surprisingly, an 
unacceptable fit (x 2 /d.o.f.=4502/1063). This re- 
sult already is an indication that variations have 
occurred in the ionized gas between the two epochs 
(2006 and 2009) and we will explore their nature 
by varying one parameter at a time in each of the 
three absorber's phases. 

• Test a): changes in column density 

We left Nfl- of the three phases free to vary 
and kept the other parameters fixed. The 
resulting model shows that the column den- 
sity of the gas in the two highest ioniza- 
tion phases (log £ ~2.3 and 3.3) go basi- 
cally to zero, but the one of the low ion- 
ization phase (log £ ~1.9) is well measured 
and it is equal to N ff =2.72±0.80x 10 20 cm" 2 
with x 2 /d.o.f=1325/1060. If the ionization 
parameter of this low ionization absorber is 
also left free, we find Iog£=1.99]to'o3 and a 
best fit of x 2 /d.o.f=1324/1059. This result 
indicates that the low-ionization phase of the 
absorber is required by the fit. The gas ion- 
ization state is consistent with the mid-state 
findings but the column density of the ab- 
sorber in 2006 is found an order of magni- 
tude lower compared to the 2009 measure- 
ment (Table E]). 

• Test b) : changes in the ionization parameter 



As a first test, we fixed the column density to 
the 2009 values and calculated the ionization 
parameters that the gas is expected to have 
if it was in photoionization equilibrium, i.e. 
we assume that the ionization state of the 
gas increases with increasing flux, and that 
the gas responds to the change of the source 
luminosity in a time scale shorter than the 
3 years elapsed between the two observa- 
tions. As reported in Table O the source 
flux changed by a factor of ^6.6 between 
the observations in mid state (2009) and in 
high state (2006). Therefore, taking into 
account the relation of the ionization pa- 
rameter to the X-ray luminosity £=L/n e r 2 
and the proportionality between luminosity 
and flux, we estimate the following ioniza- 
tion states for the three absorbers under the 
assumption of photoionization equilibrium: 
log£= 2.74, 3.17 and 4.13 (from the lowest 
to the highest ionization state, respectively). 
By forcing the gas to assume these ioniza- 
tion parameters, the fit worsens considerably 
(x 2 /d.o.f.=2518/1063). This result shows 
that a simple picture where the gas is in 
photoionization equilibrium is not consistent 
with the data. However, test a) suggests 
that some changes in column density occur. 
Therefore, to test a scenario where the gas 
is in photoionization equilibrium and at the 
same time changes in column density, we 
left Nh free in the three absorbers with ion- 
ization parameters fixed to the equilibrium 
value. None of the absorbers is robustly de- 
tected, we find upper limits of Nff<lxl0 20 , 
5xl0 21 , 2.8 xlO 20 cm" 2 from the lowest to 
the highest ionization state. However, the fit 
(x 2 /d.o.f.=1355/1060) is worse than the one 
described in case a) (where the ionization 
parameters were fixed to the 2009 values), 
indicating that at least one of the absorbers 
is not consistent with being in equilibrium. 

• Test c) : photoionization equilibrium only in 
the two highly ionized absorbers 
When the gas reaches a very high ioniza- 
tion state basically no absorption features 
are produced and the absorber is effectively 
transparent. For this reason, we cannot yet 
exclude this possibility in the two highest 
ionization components since they are always 
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undetected in the high state. To carry out 
this further test, we left the column den- 
sity and the three ionization states free to 
vary. In this way, we mainly intended to 
test the change in ionization of the low ion- 
ization component: in facts, the two higher 
phases become almost completely ionized 
and undetectable, in agreement with pho- 
toionization equilibrium, but in the low ion- 
ization absorber we found log£=1.99±0.06 
and N H =3.08±^xl0 20 cm" 2 . 

This low ionization component is the only 
one required by the fit statistics in the high 
state data. 

• Test d): changes in covering fraction 

As a final test we left the covering frac- 
tion of the three phases free to vary as it 
was done for the mid state data and de- 
scribed at the beginning of Section 13.41 The 
model used in this test is the same as the 
2009 data and the result is perhaps expected: 
the high and medium phases are basically 
transparent with covering factors equal to 
0, but the low ionization phase is required 
(x 2 /d.o.f.=1338/1057) with a covering fac- 
tor of 0.20±0.05. 

We conclude that between the high state of 
2006 and the mid state of 2009, two layers of 
the ionized gas (high and medium ionization) may 
have changed cither ionization state, in agreement 
with the prescriptions of photoionization equilib- 
rium, or they may have changed covering fraction. 
The lowest ionization phase of this gas is required 
to be present in 2009 and 2006, but when the flux 
increases it does not follow the trend of the other 
two since also in the high state data it is detected 
with approximately the same ionization state as 
of the mid state spectrum. 

3.4-3. The high state spectrum in 2000 

For completeness, we repeated the previous 
analysis on the other high state data set that is 
available, the one obtained when the source was 
observed for the first time by XMM-Newton in 
2000. This observation was considerably shorter 
compared to the one in 2006 while the flux level 
is very similar and the signal-to-noise in the RGS 
spectrum does not allow us to constrain the warm 



absorber parameters very tightly. Nonetheless, af- 
ter repeating the exercise, we can say that the 
low ionization phase of the absorber may have 
been present at this epoch with a coverage of 20%, 
therefore this spectrum is fully consistent with the 

2006 high state data. 

3.4-4- The low state spectrum in 2007 

The spectrum extracted from the observation of 

2007 corresponds to the lowest flux state observed 
so far in Mrk 335 (see the green bottom data in 
Figure [lj . For this reason, it is also the one with 
the lowest signal-to-noise in the RGS, which makes 
the analysis of warm absorber features fairly diffi- 
cult. 

We have taken the same approach described in 
the previous sections and investigate the pres- 
ence of the three absorbers in the low state 
data. At first, a model without ionized absorp- 
tion provides already an acceptable fit to the 
data, with x 2 / d -°- f -= 1477/1070. When we ap- 
plied the 3 warm absorber model with param- 
eters fixed to the mid state best fit, we find 
a worse fit (x 2 /d.o.f.=1519/1070). Then, we 
left the absorbers' column densities free to vary 
while keeping the ionization states fixed to the 
mid state values. This test provides a fit with 
X 2 /d.o.f.=1480/1067 and it shows that the col- 
umn densities of the medium and high ionization 
absorbers go to zero, while the low ionization one 
seems to survive, although it is barely detected 
with a column density of 1.67±1.10xl0 21 cm 2 . 
We calculated the upper limits of the medium and 
high ionization components and they are respec- 
tively 8xl0 20 and 1.8xl0 22 cm 2 . 

The outcome of this test shows that the low 
state data are consistent with the presence of an 
ionized absorber with a low ionization parameter, 
but this component is not formally required by the 
fit. 

It should not be surprising that a prior detailed 
analysis of the high resolution spectra of the low 
state data did not report on the presence of line- 
of sight absorption in Mrk 335 (Longinotti et al. 
2008). We recall that when the source was ob- 
served in 2007 and the data were subsequently an- 
alyzed and studied, there was no recognition of the 
presence of warm absorption in this source so that, 
given also the mcagrencss of these data, discarding 
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an in depth analysis of intrinsic absorption seemed 
a reasonable choice. 

We will discuss the implications of the results 
described in this section after presenting the re- 
sults from the HST analysis. 

4. Analysis of the HST-COS spectra 

As a first step in searching for absorption fea- 
tures in the COS spectrum that might correspond 
to the X-ray absorption observed with XMM- 
Newton, we examine the COS spectra in veloc- 
ity space as shown in Figure |U No absorption 
intrinsic to Mrk 335 is immediately obvious in 
these spectra. All narrow absorption features are 
either foreground Milky Way interstellar absorp- 
tion lines, or foreground intergalactic Lya lines. 
The foreground interstellar lines have velocities 
near the observed Milky Way Hi velocity (Mur- 
phy et al. 1996), or coincident with the Magellanic 
Stream (Fox et al. 2010). Pcnton et al. (2002) 
have identified four intergalactic Lya lines; none 
of these have counterparts in other ions. 

To carry out a more sensitive search for ab- 
sorption features in the COS spectra of Mrk 335, 
we first developed a model for the continuum and 
emission lines so that we could normalize the spec- 
trum and enhance the contrast of any intrinsic ab- 
sorption features. Since only the COS observation 
from February 2010 covered the full FUV spec- 
tral range, we started with that observation for 
our emission model. For the continuum we used a 
power law, F\ = F 10 qq (A/1000A)~ q reddened by 
E(B-V) = 0.035 (Schlegel et al. 1998) using the 
extinction law of Cardelli et al. (1989) with a ratio 
of selective to total extinction of Ry = 3.1. To get 
the best fit in the Lya and Civ regions, we op- 
timized the power-law parameters independently 
in each spectral region. We used multiple Gaus- 
sian components to trace the emission lines. For 
each line we used the minimum number of compo- 
nents to obtain a good fit. If an additional compo- 
nent improved the fit by a statistically significant 
amount (> 95% confidence using an F-test), we 
included it in our model. 

For the narrowest C iv and N v emission lines in 
Mrk 335, we include contributions from each line 
in the doublets for each component. We linked 
the doublet lines so that their wavelengths are in 
the ratio of their laboratory values. For C iv we 
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Fig. 8. — Regions in velocity space in the COS 
spectra of Mrk 335 surrounding the velocities ob- 
served for the absorption features in the X MM 
Newton mid-state spectra of Mrk 335. Relative 
fluxes are plotted as a function of velocity rela- 
tive to the systemic redshift of z = 0.025785. The 
top panel shows velocities relative to Lya in the 
COS observation of 2010 February. The middle 
panel shows the region relative to the Nv A1238 
line from the same COS observation. The bottom 
panel shows the region relative to the Civ A 1548 
line from the 2009 October COS observation. Dot- 
ted vertical lines indicate the velocities of the ab- 
sorption components fitted in the XMM-Ncwton 
spectra. 



14 



assume an optically thin ratio of 2:1 for the fluxes 
in the blue and red components of the doublets. 
The optically thin assumption did not provide a 
good fit for N v, so there we allowed the fluxes to 
vary independently. For C iv our best fit required 
an additional very broad component which is not 
needed in the fits to the Lya and Nv regions. The 
13,000 kms -1 width of this component is so much 
greater than the 500 kms -1 splitting of the Civ 
lines that we did not treat this line as a doublet. 
We used the IRAlU task specfit (Kriss 1994) to 
do our fits. 

Using this same model, we then tried to fit the 
COS G160M spectrum from October 2009. Sur- 
prisingly, although the spectra look virtually iden- 
tical at first glance, we found large residuals on the 
blue wing of the C iv emission line. As shown in 
Fig. [9l there is a shallow, but prominent, depres- 
sion from 1550-1560 A on the blue wing of the C iv 
emission line in the October 2009 COS spectrum 
that is not evident in the February 2010 COS spec- 
trum nor in the prior FOS and STIS spectra. As 
Fig. [9] shows, when all four spectra are normalized 
to the same flux level at 1500 A, the continuum 
shape and the base of the C iv line profile appear 
identical in all regards except for the depression 
at ~ 1555 A that appears in the 2009 spectrum. 
It is difficult to think of a purely emission model 
that would account for such a variable feature. We 
conclude that it is absorption intrinsic to Mrk 335. 

To rule out instrumental effects, we examined 
the four individual spectra obtained at different 
central wavelength settings on 2009 October 31 
and verified that this same feature is present in 
each one. There are a few bad pixels in this wave- 
length range in the 1611 central wavelength set- 
ting, but these are properly flagged and masked 
out in the combined spectrum. Observations of 
the white dwarf spectral standard WD1057+719 
from September and November 2009 that bracket 
the observation of Mrk 335 show only a smoothly 
varying continuum from 1540-1580 A. 

To further test our hypothesis that the spectral 
feature in the 2009 COS spectrum of Mrk 335 is 
due to absorption, we have tried several alterna- 



3 IRAF ( |http:/ /iraf.noao.edu/} is distributed by the Na- 
tional Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National 
Science Foundation. 



tive emission models for the base of the C iv emis- 
sion line and the surrounding continuum. Model 
I, our best-fit model, is the one we have described 
above. For Model II, we omit the broad 13,000 
kms -1 component of Civ that is not present in 
Lya or Nv. For Model III, we include the low- 
level forest of Fe II emission lines that are com- 
mon in AGN, particularly NLSls like Mrk 335. 
Our fit uses the relative intensities of the Wills et 
al. (1985) Fe n model, and we permit the over- 
all normalization, the line width, and the redshift 
to vary freely. In our best fit, the Fe II lines have 
FWHM=4492 ± 157 km s" 1 and a velocity relative 
to systemic of — 407±71 kms -1 . As summarized in 
Table HI Models II and III give dramatically worse 
fits. For the October 2009 COS spectrum, these 
models give A% 2 > 100, and they can be discarded 
as alternatives to our best-fit model using an ab- 
sorption feature at > 5er confidence. Again, we 
conclude that intrinsic absorption is the best de- 
scription for the unusual feature on the blue wing 
of the October 2009 COS spectrum of Mrk 335. 

For our final best-fit model (Model I) to charac- 
terize the Civ absorption in Mrk 335, we exclude 
the spectral range 1544-1565 A for our fits to the 
emission model that we apply to the COS spectra 
and the other HST spectra obtained with FOS and 
STIS. For these fits we allowed the line and contin- 
uum fluxes and line velocities and widths to vary 
freely. Our best-fit model has a power-law normal- 
ization of Fiooo = 9.43 x 10~ 14 ergs cm -2 s _1 A 
with a spectral index of a = 1.88. Table [7] gives 
the parameters of the emission-line components in 
our model. Fig. [TU] shows these best-fit emission 
models for the two COS spectra of Mrk 335 in the 
C iv region. 

For a better view of the intrinsic absorption 
in Mrk 335, we then divided our flux-calibrated 
spectra by these emission models to produce nor- 
malized spectra. Figure [TT] shows the normalized 
spectra for the Lya, N v, and C IV spectral regions 
of the two COS observations. For ease of compar- 
ison among the spectra, we have smoothed them 
with a running average of 50 kms -1 in width, or 
3 resolution elements for COS. The most promi- 
nent features in each spectrum are the foreground 
Milky Way and Magellanic Stream absorbers, as 
well as the foreground Lya IGM absorption iden- 
tified by Penton et al. (2002). In addition, in the 
October 2009 COS spectrum, we clearly see the 
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smooth, broad absorption in ClV that is intrinsic 
to Mrk 335. This feature may be weakly present in 
the February 2010 COS spectrum, but we see no 
counterparts in the Lya and Nv spectral regions. 

As confirmation that the broad, blue-shifted 
C iv absorption present in the October 2009 COS 
spectrum is real, in Figure IT21 we show normalized 
spectra of the C iv region of all HST spectra of 
Mrk 335. Again, the foreground Milky Way and 
Magellanic Stream absorbers in C iv and C I A 1560 
are the dominant features. Some weak absorption 
also appears to be present in February 2010; we 
see hints of the broad, blue-shifted C iv absorp- 
tion in the FOS spectrum, but it is absent in the 
STIS spectrum. 

To measure the column density represented by 
this absorption, or to set upper limits on its pres- 
ence, and to establish its statistical significance, 
we used two different models. Given the breadth 
of the feature, both C iv multiplcts are blended 
together, so we can not empirically integrate the 
optical depth across the spectrum to determine the 
column density and covering fraction. Instead, we 
fit a model consisting of Voigt profiles for each of 
the C iv lines, with the column density times the 
oscillator strength, Nf, set at a 2:1 ratio for the 
blue vs. the red multiplet. We tested two differ- 
ent limits — first, we fixed the covering fraction at 
unity and fit for the best position, width, and col- 
umn density; second, we set the column density to 
a highly saturated value (N = 5 x 10 19 cm" 2 ) and 
allowed the covering fraction to vary freely. To 
establish the statistical significance of the absorp- 
tion, we also performed fits with no absorption and 
allowed all emission components (line and contin- 
uum) to vary freely. Table [5] gives the best-fit pa- 
rameters for each of the HST observations. The 
blue-shifted C iv absorption is statistically signif- 
icant for both COS spectra (A%2 = -388.4 and 
-70.4) for three additional degrees of freedom for 
the October 2009 and February 2010 spectra, re- 
spectively. For the STIS and FOS spectra, the ad- 
dition of the absorption component docs not pro- 
duce a statisically significant improvement in the 
the fit, so we only quote 2a upper limits on the 
C iv column density for these spectra. 

The optically thin case gives the best fit to 
the COS spectra, which is not surprising, given 
the smooth profile of the absorption trough. The 
optically thick, partially covered model has flat- 



bottomed profiles with sharper sides. As in other 
AGN, however, it is also possible that covering 
fraction varies smoothly with velocity for an op- 
tically thick line, but characterizing a case like 
this would be highly model dependent for these 
data given the blending of the C iv multiplets. Al- 
though the optically thick case gives a significantly 
worse fit (A%2 = +31.6 for 19667 degrees of free- 
dom), it is statistically acceptable, and it allows a 
direct comparison to the X-ray models fit to the 
XMM-Newton RGS spectra that allow for partial 
covering. 

The weak detection of C iv absorption in the 
February 2010 COS spectrum suggests that cor- 
responding absorption should be present in Lya 
as well as Nv. However, each of these absorption 
troughs fall in regions of the spectrum in which it 
would be very difficult to disentangle them from 
other spectral features. Referring to Fig. [8j note 
that the Lya trough falls on the red wing of the 
damped Lya absorption from foreground Milky 
Way gas; the blended N v troughs would roughly 
be centered on the peak of the Lya emission line of 
Mrk 335. Given the breadth of the trough and its 
similar width to the broad component of the Lya 
emission line, a model for the combined absorp- 
tion and emission is highly degenerate. Our fits 
show that slight variations in the parameters for 
the Lya emission line as well as the damped Galac- 
tic Lya absorption can easily obviate the need for 
intrinsic absorption. Likewise, absorption compa- 
rable to the depth (6% of the continuum intensity) 
of the C iv absorption observed in the 2010 COS 
spectrum also provides an acceptable fit. 

For completeness, we have fit a model to the 
2010 COS spectrum that includes Lya and N V ab- 
sorption using the same velocities and widths for 
the Voigt parameters that provide the best fit to 
the October 2009 C IV absorption. Formally, the 
model that includes absorption provides a better 
fit (Ax 2 = -19.7 for 20079 degrees of freedom 
and two additional free parameters), but given 
the strong degeneracy among the components of 
the model, we simply quote 2a upper limits on 
the column densities for H I and N v of N(H i) 
< 7.7 x 10 13 cm" 2 and N(N v) < 1.0 x 10 14 cm" 2 
for the February 2010 COS spectrum. A hydro- 
gen column density as high as 1.25 x that ob- 
served for C iv in the 2010 COS spectrum, N(H i) 
= 2.1 x 10 14 cm" 2 , gives A^ 2 = +20.2 relative to 
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Table 7: Emission Features in the 2010 COS Spectra of Mrk 335 



Feature A Flux v sys FWHM 





(A) 


(10" 14 erg cm" 2 s" 1 A-1) 


(km s 1 ) 


(km s" 1 ) 


C in* 


1175.8 


18.0 ±0.5 


353 ± 55 


4727 ± 13 


Lya 


1215.6700 


16.8 ±0.6 


63 ±2 


437 ± 17 


Lya 


1215.6700 


175.0 ± 1.2 


-25 ±2 


1050 ± 15 


Lya 


1215.6700 


196.0 ±0.8 


20 ±4 


2742 ± 19 


Lya 


1215.6700 


233.0 ±0.6 


-526 ±8 


5976 ± 151 


N v 


1238.8210 


7.5 ±0.2 


-8 ±8 


1050 ± 15 


N v 


1242.8040 


7.5 ±0.2 


-9 ±8 


1050 ± 15 


N v 


1238.8210 


13.0 ±0.4 


-242 ± 4 


2742 ± 19 


N v 


1242.8040 


6.5 ±0.2 


-242 ±4 


2742 ± 19 


N v 


1238.8210 


22.6 ±0.7 


-190 ±8 


5976 ± 151 


N v 


1242.8040 


13.3 ±0.7 


-190 ±8 


5976 ± 151 


N iv] 


1486.4960 


5.5 ±0.2 


-15 ±21 


1705 ± 61 


Civ 


1548.1950 


54.9 ± 1.1 


61 ±9 


1119 ± 11 


Civ 


1550.7700 


27.5 ±0.5 


61 ±9 


1119 ± 11 


Civ 


1548.1950 


35.9 ±0.5 


-178 ±18 


2278 ± 63 


Civ 


1550.7700 


18.0 ±0.3 


-178 ±17 


2278 ± 63 


Civ 


1548.1950 


93.4 ±0.3 


-124 ±9 


5708 ± 20 


Civ 


1550.7700 


46.7 ±0.1 


-124 ±9 


5708 ± 20 


Civ 


1549.0500 


26.7 ±0.8 


-1227 ±9 


13263 ± 84 


Fe ii 


1608.4500 


30.1 ±0.4 


586 ± 71 


9992 ± 105 


He ii 


1640.4800 


8.7 ±0.2 


-103 ±7 


870 ± 19 


He ii 


1640.4800 


39.3 ±0.5 


-317 ±28 


4954 ± 74 


ml 


1663.4450 


25.6 ±0.3 


368 ± 18 


3306 ± 47 
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the best fit. Although significantly worse than the 
best fit, such a model is statistically acceptable. 
The lack of strong, obvious Lya and Nv absorption 
features in the 2010 COS spectrum is not unusual 
given the weakness of the C iv absorption at this 
epoch. 

5. Discussion 

The most recent X-ray and UV observations of 
Mrk 335 here presented have revealed an unex- 
pected behavior of this complex source. 

These observations have provided us with the 
opportunity to study the emergence of an out- 
flow in a source where this phenomenon was not 
observed. Among similar cases reported in the 
literature, we recall the two luminous quasars 
J105400.40+034801.2 (Hamann et al. 2008) and 
Ton 34 (Krongold et al. 2010), and the Narrow 
Line Scyfert 1 WVPS 007 (Leighly et al. 2009 and 
reference therein) . In all cases the outflow was re- 
vealed in the UV broad absorption lines, a fairly 
common phenomenon in luminous quasars. The 
case of WVPS 007 is different since broad absorp- 
tion is not expected in low luminosity objects such 
as this Seyfert 1 type AGN (see Laor & Brandt, 
2002). Furthermore, in WVPS 007 the UV out- 
flow is accompanied by extreme X-ray variability 
(Grupe et al. 2008b). The new findings on Mrk 
335 presented in this work arc reminiscent of the 
behavior of WVPS 007, with the advantage of of- 
fering the "X-ray view of the absorber" thanks to 
its higher X-ray flux. 

Warm absorber variability has been extensively 
observed and studied in several Seyfert 1 Galaxies, 
but in most of the reported cases, the ionized gas 
has been a constant characteristic of the source, 
albeit with changing properties. 

The presence of ionized gas in Mrk 335 is some- 
what elusive in its X-ray records. The results from 
the ASCA survey of AGN published by Reynolds 
(1997) excluded the presence of ionized absorption 
in Mrk 335 spectra. Subsequently, observations 
by high throughput facilities tend all to confirm 
the "bareness" of the intrinsic line of sight to this 
Seyfert Galaxy (Bianchi et al. 2001, O' Neill et al. 
2007, Larsson et al. 2007). 

Quite remarkably, the presence of a variable 
ionized absorber uncorrelated with the source flux 
was first proposed in the earlier decade to explain 



the Ginga spectra of Mrk 335 (Turner et al. 1993). 
These authors detected a column density around 
1.4xl0 22 in the Ginga spectrum of 1987, which 
was inconsistent with the Ginga data taken the 
following year. 

Since then, we had to wait until the XMM- 
Ncwton visit of 2009 to observe the ionized gas 
crossing our line of sight again. The availabil- 
ity of multiple XMM-Newton and Swift spectra 
allowed Grupe and collaborators (2012, and ref- 
erence therein) to propose partial covering from 
intervening gas as a viable explanation for most of 
the strong variability exhibited by Mrk 335, albeit 
an alternative interpretation in terms of blurred 
reflection remained equally valid to explain the 
CCD spectra (see also Gallo et al. 2013). Luckily, 
this time we could rely on the observables mea- 
sured by a high resolution spectrometer and re- 
veal the properties of the ionized gas with unprece- 
dented detail. 

The analysis presented here favors a solution 
where the variability of the absorber is related 
to a change in the column density and/or in the 
covering fraction of the gas. Changes of the ab- 
sorber's properties correlated with the flux vari- 
ability have to be discarded on the basis of Test 
b) in section [3.4.21 This test has shown that the 
ionization parameter expected in photoionization 
equilibrium for the low ionization component of 
the absorber is not consistent with the data. We 
avoided carrying out the same test in the low state 
data (scction|3A4|) since it gave already a negative 
response in the data set with the highest signal to 
noise: if the ionization parameter was responding 
to the flux increase, this change would be more 
likely detected in the RGS spectra of 2006 and 
our results in sect ion [3 . 4 . 2 1 have shown that this is 
not the case. 

This lack of correlation therefore points out 
that the warm absorber has overall little effect on 
the variability of the continuum. 

The results from the HST analysis provide us 
with additional information on the intriguing be- 
havior of this source. Our photoionization models 
for the X-ray warm absorber in Mrk 335 predict 
substantial columns of UV-absorbing gas (column 
densities for Lya, Nv, and Civ are respectively 
200, 100 and 30 xlO 14 cnT 2 for the low ionization 
component of the X-ray absorber). The column 
densities measured using COS (see Table |8]) are 



19 



Table 8: Best Fit Parameters for C IV Absorption in the HST Spectra of Mrk 335 



Instrument 


Date 


Model 3 


N(Civ) 


v sys 


FWHM 


fcov 


x'/d.o.f. 




(cm -2 ) 


(km s" 1 ) 


(km s _1 ) 










FOS 


1994-12-16 


I 


< 0.8 x 10 14 


-6503 (fixed) 


2602 (fixed) 


1.0 (fixed) 


875.26/778 


FOS 


1994-12-16 


1 


5 x 10 19 (fixed) 


-6503 (fixed) 


747 (fixed) 


0.03 ±0.01 


876.08/778 


FOS 


1994-12-16 


1 


0.0 (fixed) 








881.52/779 


STIS 


2004-07-01 


1 


< 0.5 x 10 14 


-6503 (fixed) 


2602 (fixed) 


1.0 (fixed) 


19696.46/18666 


STIS 


2004-07-01 


1 


5 x 10 19 (fixed) 


-6503 (fixed) 


747 (fixed) 


< 0.02 


19695.35/18666 


STIS 


2004-07-01 


1 


0.0 (fixed) 








19676.12/18667 


COS 


2009-10-31 


1 


(4.9 ±0.1) x 10 14 


-6503 ± 34 


2602 ± 57 


1.0 (fixed) 


13756.76/19668 


COS 


2009-10-31 


I 


5 x 10 19 (fixed) 


-6503 ± 14 


747 ±7 


0.15 ±0.003 


13788.39/19669 


COS 


2009-10-31 


I 


0.0 (fixed) 








14368.39/19671 


COS 


2009-10-31 


11 


0.0 (fixed) 








14313.10/19674 


COS 


2009-10-31 


111 


0.0 (fixed) 








14467.21/19671 


COS 


2010-02-08 


1 


(1.7 ±0.1) x 10 14 


-6010 ±49 


1744 ± 181 


1.0 (fixed) 


13565.24/19668 


COS 


2010-02-08 


1 


5 x 10 19 (fixed) 


-6494 ± 38 


731 ± 25 


0.06 ± 0.004 


13601.12/19667 


COS 


2010-02-08 


1 


0.0 (fixed) 








13635.63/19671 


COS 


2010-02-08 


11 


0.0 (fixed) 








13621.51/19674 


COS 


2010-02-08 


III 


0.0 (fixed) 








13808.96/19671 



a Model I uses a power-law continuum and a broad base on the C IV emission line, and C IV absorption is 
treated as a Gaussian in optical depth. Model II uses a power-law continuum and no broad base for the C IV 
emission line. Model III uses a power-law continuum and a Wills et al. (1985) model for Fe II emission. Best 
fits for this model prefer no broad base on the C IV emission line. 



comparable to these predicted column densities, 
suggesting that the UV absorption we see may 
be related to the X-ray absorption. The veloci- 
ties of the UV absorption troughs are also roughly 
comparable to the X-ray absorption, but we note 
that the UV absorption extends to higher veloci- 
ties and covers a larger total range. In this exercise 
of comparing the X-ray and UV absorption, we al- 
ways need to keep in mind that the observations 
are non simultaneous, hence the effect of intrinsic 
variations of these column densities that may oc- 
cur in between June and (late) October 2009 and 
February 2010, cannot be taken into account. 

Our COS observations show that substantial 
variations in the C IV absorption can occur on 
time scales of a few months as well as on time 
scales of years. If the change in absorption is due 
to a change in column density, this could be due 
to either bulk motion of gas or due to an ioniza- 
tion response. However, between October 2009 
and February 2010, there is very little change in 
the UV continuum flux. Likewise, longer-term UV 
continuum monitoring using Swift shows that UV 
flux variations in Mrk 335 are only at the level of 
tens of percent (Grupe et al. 2012). These contin- 
uum variations are not strong enough to generate 
column-density variations as strong as we see. 




Wavelength (A) 



Fig. 10.— The two COS spectra of Mrk 335 with 
best-fit models of the continuum and line emis- 
sion overlayed. The 2009 October 31 spectrum is 
in black, and the 2010 February 8 spectrum is in 
red. We have smoothed the spectra with a running 
average of 50 km s _1 in width. The blue and green 
lines show the fitted power-law continuum and the 
total fit including all emission-line components. 
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Fig. 11. — Normalized spectra of Mrk 335 in veloc- 
ity space. Normalized relative fluxes are plotted 
as a function of velocity relative to the systemic 
rcdshift of z = 0.025785. The spectra have been 
smoothed with a running average of 50 kms _1 in 
width. Dotted vertical lines indicate the veloci- 
ties of the absorption components in the XMM- 
Newton mid-state spectra. The top panel shows 
the Lya region of the COS spectrum of 2010-02- 
08. Here, the IGM Lya lines identified by Pen- 
ton et al. (2002) are marked. The second panel 
shows the Nv region of the 2010-02-08 spectrum 
with Nv A1238 in blue and Nv A1242 line in red. 
The two bottom panels show the C iv region of the 
COS spectra of 2010-02-08 and 2009-10-31, respec- 
tively; Civ A1548 is in blue, and A1550 is in red. 
The narrow absorption features in each panel are 
foreground Galactic and IGM lines, which we have 
labeled. 



-10000 -5000 

Velocity (km s _1 ) 

Fig. 12. — Same as Figure [11] but including all 
archival HST data. From top to bottom: FOS 
spectrum of 1994-12-16, STIS spectrum of 2004- 
07-01, COS spectrum of 2009-10-31 and COS spec- 
trum of 2010-02-08. To facilitate the comparison 
of the spectra from instruments of differing reso- 
lution, we have smoothed the spectra with a run- 
ning average of 50 km s _1 in width. Dotted vertical 
lines indicate the velocities of the absorption com- 
ponents in the XMM-Newton mid-state spectra. 
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Therefore, we favor bulk motion of the gas as 
the best explanation for the absorption variations, 
which is nicely in agreement with the results from 
the X-ray analysis. 

If the X-ray and UV absorbing gas are rapidly 
variable and ephemeral, as in Mrk 766 (Risaliti et 
al. 2011), the absorbers may have moved trans- 
versely across our line of sight in the interven- 
ing time between the COS observations. Assum- 
ing that such motions are due to Keplerian orbits 
around the central black hole in Mrk 335, we can 
set limits on the distance of the outflowing gas 
from the central black hole. For a black hole mass 
of 2.6 x 10 7 Mq (Grier et al. 2012), the minimum 
size for a cloud that fully covers a disk emitting 
region of 100 Rg is 4.2 x 10 14 cm. For it to move 
across our line of sight in the 100 days between the 
two COS observations, the velocity would have to 
be only 450 kms -1 , which would place it at a dis- 
tance of 3.2 x 10 18 cm, slightly exterior to the torus 
and in the inner narrow-line region (NLR). But, at 
that distance, it would also cover the broad-line re- 
gion (BLR), which is a few thousand Rq in size. 
To move across the BLR in that time, the trans- 
verse motion would have to be 10-100 times more 
rapid, more consistent with the outflow velocities 
of 5000 - 9000 kms" 1 seen for the X-ray and UV 
absorbers. A velocity of 5000 kms -1 would be 
in Keplerian motion at 1.4 x 10 16 cm. The BLR 
in Mrk 335 is at 0.7 - 4.7 x 10 16 cm (Grier et al. 
2012, Peterson et al. 2004), so the absorbers would 
then be in the inner BLR, as suggested for Mrk 766 
(Risaliti et al. 2011). 

This region was also proposed as a likely lo- 
cation for the X-ray emission lines seen in the 
low-state spectrum of Mrk 335 (Longinotti et al. 
2008). The new data in the mid state show that 
the fluxes of the emission lines do not vary sig- 
nificantly (see section [3~3]) . Notably, in 2009 we 
see again a dominant intcrcombination compo- 
nent with respect to the forbidden line, which 
seems to confirm the high density of the emit- 
ting gas previously postulated (Longinotti et al. 
2008). Whether the emitter and the absorber 
are physically related is difficult to assess without 
more stringent constraints on the absorber loca- 
tion. The limits on the column densities of the 
emitting gas derived by these authors in the low 
state data are large enough to be consistent with 
the values measured in this work for the absorbers. 



High- velocity transverse motion is a natural ex- 
planation for the differences that we have seen in 
the UV spectra. UV absorbers originating interior 
to the BLR would be much closer to the central 
source than has been found to date for most UV 
absorbers, which are more commonly found in the 
inner narrow-line region (Crenshaw et al. 2005). 
Such clouds with transverse and outflow velocities 
of thousands of kms -1 would be the first indica- 
tion of material at the base of a disk wind. A cloud 
moving transversely at 5000 kms" 1 would cross 
the radiating accretion disk in 10 days. Therefore, 
we would expect that future monitoring observa- 
tions should see significant absorption variations 
on timescales of weeks. As a final remark, we high- 
light once more that deep simultaneous UV and 
X-ray data will be beneficial for testing the com- 
pelling results presented in this paper and charac- 
terizing the outflow with higher detail. 
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